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The quantum-limited linewidth of a laser cavity is enhanced above the Schawlow-Townes value by 
the Petermann factor K, due to the non- orthogonality of the cavity modes. We derive the relation 
between the Petermann factor and the residues of poles of the scattering matrix and investigate 
the statistical properties of the Petermann factor for cavities in which the radiation is scattered 
chaotically. For a single scattering channel we determine the complete probability distribution of 
K and find that the average Petermann factor {K) depends non-analytically on the area of the 
opening, and greatly exceeds the most probable value. For an arbitrary number A'' of scattering 
channels we calculate (A") as a function of the decay rate F of the lasing mode. We find for ^ 1 
that for typical values of F the average Petermann factor {K) oc \fN 3> 1 is parametrically larger 
than unity. 

PACS; 42.50.Lc, 42.5Q.Ar, 42.60.Da 
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I. INTRODUCTION 

Laser action selects a mode in a cavity and enhances 
the output intensity in this mode by a non-hnear feedback 
mechanism. Vacuum fluctuations of the electromagnetic 
field ultimately limit the narrowing of the emission spec- 
trum The quantum-limited linewidth, or Schawlow- 
Townes linewidth. 
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(1.1) 



is proportional to the square of the decay rate F of the 
lasing cavity mode and inversely proportional to the out- 
put power / (in units of photons/s). This is a lower 
bound for the linewidth when F is much less than the 
linewidth of the atomic transition and when the lower 
level of the transition is unoccupied. Many years after 
the work of Schawlow and Townes it was realized 
that the true fundamental limit is larger than Eq. 
by a factor K that characterizes the non-orthogonality of 
the cavity modes. This excess noise factor, or Petermann 
factor, has generated an extensive literature 

Apart from its importance for cavity lasers, the Peter- 
mann factor is of fundamental significance in the more 
general context of scattering theory. A lasing cavity 
mode is associated with a pole of the scattering matrix 
in the complex frequency plane. We will show that the 
Petermann factor is proportional to the squared modulus 
of the residue of this pole. Poles of the scattering ma- 
trix also determine the position and height of resonances 
of nuclei, atoms, and molecules [^. Powerful numerical 
tools that give access to poles even deep in the complex 
plane have been developed recently They can be 

used to determine the residues of the poles as well. Our 
work is of relevance for these more general studies, be- 
yond the original application to cavity lasers. 




FIG. 1. Chaotic cavity that radiates light from a small 
opening. 



Existing theories of the Petermann factor deal with 
cavities in which the scattering is essentially one- 
dimensional, because the geometry has a high degree of 
symmetry. For such cavities the framework of ray optics 
provides a simple way to solve the problem in a good ap- 
proximation 1^. This approach breaks down if the light 
propagation in the cavity becomes chaotic, either because 
of an irregular shape of the boundaries (like for the cavity 
depicted in Fig. |^) or because of randomly placed scatter- 
ers. The method of random-matrix theory is well-suited 
for such chaotic cavities Instead of considering 

a single cavity, one studies an ensemble of cavities with 
small variations in shape and size, or position of the scat- 
terers. The distribution of the scattering matrix in this 
ensemble is known. Recent work has provided a detailed 
knowledge on the statistics of the poles [T^-[l^. Much 
less is known about the residues ||r^,|l^ . In this work we 
fill the remaining gap to a considerable extent. 

The outline of this paper is as follows. In Section || 
we derive the connection between the Petermann factor 
and the residue of the pole of the lasing mode. The 
residue in turn is seen to be characteristic for the de- 
gree of nonorthogonality of the modes. In this way we 
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make contact with the existing literature on the Peter- 
mann factor ||7|,p| . 

In Section we study the single-channel case of a 
scalar scattering matrix. This applies to a cavity that 
is coupled to the outside via a small opening of area 
A ^ /2tt (with A the wavelength of the lasing mode). 
For preserved time-reversal symmetry (the relevant case 
in optics) we find that the ensemble average oi K — 1 
depends non- analytically oc TlnT~^ on the transmission 
probability T through the opening, so that it is beyond 
the reach of perturbation theory even if T ^ 1. We 
present a complete resummation of the perturbation se- 
ries that overcomes this obstacle. We derive the condi- 
tional distribution P{K) of the Petermann factor at a 
given decay rate F of the lasing mode, valid for any value 
of T. The most probable value oi K — 1 is cx T, hence it 
is parametrically smaller than the average. 

In a cavity with such a small opening the deviations 
of K from unity are very small. For larger deviations we 
study, in Section IV , the multi-channel case of an x 



scattering matrix, which corresponds to an opening of 
area A « N}? /2tt. The lasing mode acquires a decay 
rate F of order Fq = NTA/2tt (with A the mean spacing 
of the cavity modes). We compute the mean Petermann 
factor as a function of F for broken time-reversal sym- 
metry, which is technically simpler than the case of pre- 
served time-reversal symmetry, but qualitatively similar. 
We find a parametrically large mean Petermann factor 
K (X Vn. 

Our conclusions are given in Section 0. The main 
results of Sections III and IV have been reported in Refs. 
|l9| and pffl, respectively. 



A small opening in the cavity is described by a real, 
non-random M x N coupling matrix W, with N the 
number of scattering channels transmitted through the 
opening. (For an opening of area yi, ~ 2ttA/X'^ at 
wavelength A.) Modes of the open cavity are complex 
eigenvalues (with negative imaginary part) of the non- 
Hermitian matrix 



n = H 



tWWK 



(2.3) 



In absence of amplification or absorption, the scatter- 
ing matrix S at frequency to is related to H. by P,El| 



s = t-2mw\uj -ny^w. 



(2.4) 



The scattering matrix is a unitary (and symmetric, for 
13 = 1) random N x N matrix, with poles at the eigen- 
values oiH. It enters the input-output relation 



N 

^(^) = ^5™„(^)aj?(t^), 



(2.5) 



which relates the annihilation operators a°"* of the scat- 
tering states that leave the cavity to the annihilation op- 
erators ajf of states that enter the cavity. The indices n, 
m label the scattering channels. 

We now assume that the cavity is filled with a homo- 
geneous amplifying medium (constant amplification rate 
I/tq over a large frequency window Qa — LA, L N). 
This adds a term i/2Ta to the eigenvalues, shifting them 
upwards towards the real axis. The scattering matrix 



S = l- 2mW\uj -H- i/2Tar^W 



(2.6) 



II. RELATIONSHIP BETWEEN PETERMANN 
FACTOR AND RESIDUE 



is then no longer unitary, and the input-output relation 
changes to p2|,p3[ 



Modes of a closed cavity, in the absence of absorp- 
tion or amplification, are eigenvalues u;„ of a Hermitian 
operator H. This operator can be chosen real if the sys- 
tem possesses time-reversal symmetry (symmetry index 
/3 = 1), otherwise it is complex (/? = 2). For a chaotic 
cavity, H can be modeled by an M x M Hermitian matrix 
with independent Gaussian distributed elements. 
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(For (3 = 1 (2), this is the Gaussian orthogonal (unitary) 
ensemble [Q.) The mean density of eigenvalues is the 
Wigner semicircle 



p{uj) 



M 
27r/x^ 



(2.2) 



The mean mode spacing at the center w = is A = 
■Kfi/M . (The limit M — > cx) at fixed spacing A of the 
modes is taken at the end of the calculation.) 



N 

\u;) = J2 Sran{u^)a':{iu) + ^ Q™„ (a;)6t (a;) . (2.7) 



ri=l 



n=l 



All operators fulfill the canonical bosonic commutation 
relations [a„(ijj), aj,j(a;')] = SnmS{uj — lj'). As a conse- 
quence. 
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(2.8) 



The operators b describe the spontaneous emission of 
photons in the cavity and have expectation value 



{bl{uj)bmil^')) = SnmS{uJ - Uj')f{uJ, T), 



(2.9) 



with f{uj,T) = [exp(?icj/fcBT) — 1] the Bose-Einstein 
distribution function at frequency lj and temperature T. 



In the absence of external illumination ((a 
the photon current per frequency interval. 
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(2.10) 
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is related to the scattering matrix by Kirchhoff's law 
11,11 

= /(a>,r)-Ltr [1 - SHu;)S{u;)]. (2.11) 

ZTT 

For Lo near the laser transition we may replace / by the 
population inversion factor N^p/ (Mow — -^up), where N^p 
and -/Viow are the mean occupation numbers of the upper 
and lower levels of the transition. In this way the photon 
current can be written in the form 
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low 



-ti[S^uj)S{uj)~l], (2.12) 



that is suitable for an amplifying medium. (Alternatively, 
one can associate a negative temperature to an amplify- 
ing medium.) 

The lasing mode is the eigenvalue Q — ir/2 closest to 
the real axis, and the laser threshold is reached when the 
decay rate T of this mode equals the amplification rate 
I/tq. Near the laser threshold we need to retain only the 
contribution from the las ing mode (say mode number I) 
to the scattering matrix (^.6|), 
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{W^UU{U-^W)i„ 



LO - 



VL + iV/2-il2Ta 



(2.13) 



where U is the matrix of right eigenvectors of TL (no 
summation over / is implied). The photon current near 
threshold takes the form 



2tiN, 



up 



{Www'^u)u{u-^ww'^u-^'^)i 



(2.14) 



This is a Lorentzian with full width at half maximum 

6uj - l/Ta- 

The coupling matrix W can be eliminated by writing 



- Tiiu^ww'^uYi = im {u^nu)^^ 

-n{U-^WW^U-^^)u =Im ([/-^WC/^^t) 



(2.15a) 



= -,^{U-'U-'%. (2.15b) 

The total output current is found by integrating over 
frequency, 



{U^U)u{U'^U~^^)i 



Mow SUJ ' 



(2.16) 



Comparison with the Schawlow-Townes value (1.1 ) shows 
that 



6uj = 2K- 



Nup 



-5uj. 
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Nup — Mow 

where the Petermann factor K is identified as 



(2.17) 



K = iU^U)uiU-^U-'%>l. (2.18) 

For time-reversal symmetry, we can choose = U'^ , 
and find K = [(f7[/^)n]^. The factor of 2 in the relation 
between Suj and Sujgj, occurs because we have computed 
the laser linewidth in the linear regime just below the 
threshold, instead of far above the threshold. The ef- 
fect of the non-linearities above threshold is to suppress 
the amplitude fluctuations while leaving the phase fluc- 
tuations intact ||2^, hence the simple factor of two re- 
duction of the linewidth. The factor A^up/(-^up — Mow) 
accounts for the extra noise due to an incomplete pop- 
ulation inversion. The remaining factor K is due to the 
non-orthogonality of the cavity modes |^ , since K ^ 1 if 
U is unitary. 



III. SINGLE SCATTERING CHANNEL 



The relation ( 2.18 ) serves as the starting point for 
a calculation of the statistics of the Petermann factor 
in an ensemble of chaotic cavities. In this section we 
consider the case N = 1 of a single scattering chan- 
nel, for which the coupling matrix W reduces to a vec- 
tor d = (IVii, W21, ■ ■ • , Wmi). The magnitude |ap = 
{MA/tt'^)w, where w € [0,1] is related to the trans- 
mission probability T of the single scattering channel by 
T = Aw{l + w)~^ We assume a basis in which H is 
diagonal (eigenvalues w^, right eigenvectors \q), left eigen- 
vectors {q\). In this basis the entries aq remain real for 
/3 = 1, but become complex numbers for (3 = 2. Since the 
eigenvectors |g) point into random directions, and since 
the fixed length of a becomes an irrelevant constraint in 
the limit M — > c», each real degree of freedom in aq is 
an independent Gaussian distributed number fl^ . The 
squared modulus |agp has probability density 
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2^|«,P 



2Ac^q? 

wA 



P/2 



exp 



2u;A' 



(3.1) 



Eq. (3.1) is a x^-distribution with f3 degrees of freedom 
and mean Aw/tt^. 

We first determine the distribution of the decay rate F 
of the lasing mode, following Ref. [Q. Since the lasing 
mode is the mode closest to the real axis, its decay rate 
is much smaller than the typical decay rate of a mode, 
which is ~ TA. Then we calculate the conditional dis- 
tribution and mean of the Petermann factor for given F. 
The unconditional distribution of the Petermann factor 
is found by folding the conditional distribution with the 
distribution of F, but will not be considered here. 



A. Decay rate of the lasing mode 

The amplification with rate l/Ta is assumed to be ef- 
fective over a window fla = LA containing many modes. 
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The lasing mode is the mode within this window that 
has the smallest decay rate F. For such small decay rates 
we can use first-order perturbation theory to obtain the 
decay rate of mode q, 



Tq = 27r|aq| 



(3.2) 



The x^-distribution (3^) of the squared moduli \aq 



translates into a x^-distribution of the decay rates, 



P(F) oc f(2-«/2 



exp 



Aw A J 



(3.3) 



Ignoring correlations, we may obtain the decay rate 
of the lasing mode by considering the L decay rates as 
independent random variables drawn from the distribu- 
tion P{T). The distribution of the smallest among the L 
decay rates is then given by 



1- [ dF'P(F') 

"'0 



L-l 



(3.4) 



For small rates F we can insert the distribution (3.2) and 
obtain 



L-l 



, = 
(3.5a) 
(3.5b) 



Here erf(a;) = 2tt^^^^ dy exp{—y^) is the error func- 
tion. The decay rate of the lasing mode decreases with 
increasing width of the amplification window as F ~ 
wA{na/A)-y^ < wA. 



B. First-order perturbation theory 

If the opening is much smaller than a wavelength, then 
a perturbation theory in a seems a natural starting point. 
We assign the index / to the lasing mode, and write the 
perturbed right eigenfunction \l)' — J2q^q\l) ^^'^ 
perturbed left eigenfunction — ^^^^^(gl, in terms of 
the eigenfunctions of H. The coefficients are dq = Uqi/Uu 
and Cq = Ur^/UiJ^, i.e., we do not normalize the per- 
turbed eigenfunctions but rather choose di = ei = 1. 

To leading order the lasing mode remains at Q — uji 
and has width 



F = 27r|az|^. 
The coefficients of the wave function are 



naaa 



na ai 
I . 

Wq - UJl 



(3.6) 



(3.7) 



The Pete rmann factor of the lasing mode follows from 
Eq. (ra, 



K 



{^+Eq^l\dq\') {l+Eq^l\eq 



1 + E 



q^l dqCq 



q^l 



(3.8) 



where we linearized with respect to F because the lasing 
mode is close to the real axis. From Eq. (3.7) one finds 



q^l 



{UJI - UJqY ' 



(3.9) 



We seek the distribution P{K) and the average {K)^ y 
of K for a given value of Vl and F. 

For /3 = 1, the probability to find an eigenvalue at 
given that there is an eigenvalue at uji vanishes linearly for 
small \ujq — LOi\, as & consequence of eigenvalue repulsion 
const rain ed by time-reversal symmetry. Since the expres- 
sion ( |3.9[ ) for K diverges quadratically for small \ujq — oj/l, 
we conclude that {K)^^y does not exist in perturbation 
theory]^ This severely complicates the problem. 



C. Summation of the perturbation series 

To obtain a finite answer for the average Petermann 
factor we need to go beyond perturbation theory. By a 
complete summation of the perturbation series we will 
in this section obtain results that are valid for all values 
T < 1 of the transmission probability. Our starting point 
are the exact relations 



dqZl — UJqdq — ilTUq ^ a*dp, 

p 



(3.10a) 
(3.10b) 



between the complex eigenvalues Zq of TL and the real 
eigenvalues LOq oi H. Distinguishing between q — I and 
q ^ I, we obtain three recursion relations. 



id„ 



7ra„ 



Zl - LOq 

TTQ!* 
Zl - UJq 



ai 



(3.11a) 
(3.11b) 
(3.11c) 



* For broken time-reversal symmetry there is no diver- 
gence. We can use the known two-point correlation func- 
tion R{u)i,ujq) of the Gaussian unitary ensemble to obtain 
{K)ns = 1 + IttTF/A for T < 1. 
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We now use the fact that zi is the eigenvalue closest to 
the real axis. We may therefore assume that zi is close 
to the unpertur bed va. lue w/ and replace the denominator 
Zi — ujq in Eq. ( 3.11c ) by loi — ujq. That decouples the 
recursion relations, which may then be solved in closed 
form, 



zi = ui — zTrjaip (1 + ittA) 



idq = 



ie„ 



LUl - LOq 



(l + iTryl) \ 



LOl — UJq 

We have defined 

q^l 

The decay rate of the lasing mode is 

r = -2Imz, ^ 2Ti\ai\^ {I + TT^ A^y^ . 



From Eq. (3.8) we find 



K=l 



2-kT 



B 



A l + 7r2v42' 



with 



B^/\Y,\aq\'{ui-UJq)- 
q^l 



(3.12a) 
(3.12b) 

(3.12c) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



The problem is now reduced to a calculation of the 
joint probability distribution P{A,B). This problem is 
closely related to the level curvature problem of random- 



matrix theory 
Appendix ^. T 

PiA,B) = 



The calculation is presented in 
le result is 



{n^A^+wY 



24 V TTW y B2+3/3/2 

2B 



X exp 



(3.17) 



D. Probability distribution of the Petermann factor 



From Eqs. (^, (gJJ), ( |3.15D , and ( |3.17D we can com- 
pute the probability distribution 



K - 1 



p{K) = {zy' (^5 
z = s{n-LJi)s I r- 



27rr B 
^ "A~ 1 + 7r2A2 

27r|a,|2 



1 + 7r2A2 / ' 



(3.18a) 
(3.18b) 



of K at fixed F and by averaging over \ai\'^, A, and B. 
In principle one should also require that the decay rates 



0, 




FIG. 2. Probability distribution of the rescaled Petermann 
factor K = (JC- 1) A/FT for T = 1 and T < 1, in the presence 
of tim e-reversal symmet ry. Th e solid curves follow from Eqs. 
( |3.20| ) (with /? = 1) and ( |3.21a[ ). The data points follow from a 
numerical simulation of the random-matrix model. The inset 



shows the results (3.20) (with /3 = 2) and (3.21b) for broken 
time-reversal symmetry. 



of modes q ^ I are bigger than F, but this extra condition 
becomes irrelevant for F — > 0. The average of Z over |a; p 
with Eq. (|^) yields a factor (1 7r2A2)^/2. (Only the 
behavior of P(|q!;|2) for small jajp matters, because we 
concentrate on the lasing mode.) After integration over 
B the distribution can be expressed as a ratio of integrals 
over A, 



P{K) 



(27r)2/3 Aw / (j^- 1)A 

" (1^^^A2/wY 
(H-^2^2)l+/3 

(l+7r2A2)'3/2 



dA 



(l+7r2A2/u)2)l+/3/2 



-2-3/3/2 



/37rwF(l +7r2A2/u;2) 

'(i^- l)A(H-7r2A2) 

-1 

(3.19) 



We introduce the rescaled Petermann factor k = {K — 
1)A/Fr. A simple result for F(k) follows for T = 1, 



P{k) 

and for T < 1, 



4/37r' 



2/3 



3^2+3/3/2 



exp 



TT 

2k 



exp 



P(k) 



/2k5 



27r 
3k 



TT 

3^ 





TT 








TT ■ 












4^. 








TT " 


1 exp 


-~2^. 



(3.20) 

(3.21a) 

/3 = 2. 
(3.21b) 



As shown in Fig. |2| the distributions are very broad and 
asymmetric, with a long tail towards large k. 
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T r/To 



FIG. 3. Average of the rescaled Petermann factor k as a 
function of tran smission probability T. The sohd curve is 
the result (3.22) in the presen ce of time-reversal symmetry, 
the dashed curve is the result (3.24) for broken time-reversal 
symmetry. For small T, the solid curve diverges oc InT"^ 
while the dashed curve has the finite limit of 7r/3. For T = 1 
both curves reach the value 27r/3. 



To check our analytical results we have also done a 
numerical simulation of the random-matrix model, gen- 
erating a large numb er of random matrices H and com- 
puting K from Eq. ( 2.18 ). As one can see from Fig. ||, 
the agreement with the theoretical predictions is flawless. 



FIG. 4. Decay-rate distribution P{V) of a chaotic cavity 
with an opening that supports A'' = 2, 4, 6, 8, 10, 12 fully 
transmitted scattering channels. Computed from Eq. (4.2), 
for the case of broken time-reversal symmetry. 



symmetry. If time-reversal symmetry is broken, then the 
stronger quadratic repulsion is sufficient to overcome the 
Lo~'^ divergence of perturbation theory (3.9) and the av- 
erage K becomes an an alytic function of T. For this case, 
we find from Eq. (3.19) the mean Petermann factor 



{K)n.r = 1 + 
shown dashed in Fig. ||. 



r 4:TTW 

A3(l-hw2)' 



(3.24) 



E. Mean Petermann factor 



The distribution (3.19) gives for preserved time- 
reversal symmetry (/? = 1) the mean Petermann factor 



(K) 



r 27r 





_ 1 _ 1 

2 2 



A 3 



^22 



1 1 

' 2 2 

-1 



(3.22) 



in terms of the ratio of two Meijer G-functions. We have 
plotted the result in Fig. 0, as a function of T = Aw{l + 
w)-\ 

It is remarkable that the average K depends non- 
analytically on T, and hence on the area of the opening. 
(The transmission probability T is related to the area A 
of the opening by T ~ A^/X^ for T < 1 H].) For T < 1, 
the average approaches the form 



(K) 



ttTT , 16 



(3.23) 



The most probable (or modal) value of if — 1 ~ TV / A 



is parametrically smaller than the mean value (3.23) for 
T <C 1. The non-analyticity results from the relatively 
weak eigenvalue repulsion in the presence of time- reversal 



IV. MANY SCATTERING CHANNELS 

For arbitrary number of scattering channels N the cou- 
pling matrix is an M x iV rectangular matrix. The 
square matrix ttW^W has N eigenvalues (MA/7r)w„. 
The transmission coefficients of the eigenchannels are 



4Wn 



(4.1) 



A single hole of area ^ ^ (at wavelength A) cor- 
responds to ~ 27r^/A^ fully transmitted scattering 
channels, with all T„ = Wn = 1 the same. 

As in the single-channel case, we first determine the 
distribution of the decay rate F of the lasing mode. This 
decay rate is smaller than the typical decay rate Fq — 
TNA/27T of the non-lasing modes. Then we calculate the 
mean Petermann factor (K) for given F and investigate 
its behavior for the atypically small decay rates of the 
lasing mode. 
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A. Decay rate of the lasing mode 

The distribution of decay rates P{T) has been cal- 
culated by Fyodorov and Sommers. For broken time- 
reversal symmetry the result is 



(4.2a) 





(- 












a 




OQ 








dx 







N 



N 



(4.2b) 
(4.2c) 



where g„ = — 1 + 2/T„. For identical gn = g the two 
functions J-i and J-2 simplify to 



My) = 



1 



(iV-l)! 

N 



(4.3a) 



n=0 

and a convenient form of the distribution function is 
2^r2(7V^/v(i_^)r/ro 



p(r) = 



(4.4) 



The behavior of P(r) for various numbers N of fully 
transmitted (T = 1) scattering channels is illustrated in 
Fig. |. 

The result for preserved time-reversal symmetry is a 
bit more involved [0. Fortunately, we can draw all 
important conclusions from the results for broken time- 
reversal symmetry, on which we will concentrate here. 

For large N, the distribution P(r) becomes non-zero 
only in the interval Fq < F < ro/(l — T), where it is 
equal to pC|,|3l|l 



P(F) 



Fo <r<ro/(i-T). 



(4.5) 



This limit is /3-independent. The smallest decay rate Fq 
corresponds to the inverse mean dwell time in the cavity. 

We are interested in the "good cavity" regime, where 
the typical decay rate Fq is small compared to the ampli- 
fication bandwidth fla- From Fq = TNA/2tt it follows 
that the number L ~ f^a/A of amplified modes is then 
much larger than TN. In this regime the decay rate of the 
lasing mode (the smallest among the L decay rates in the 
frequency window Qa) drops below Fq. The asymptotic 
result ( |4.5| ) cannot be used in this case, since it does not 
describe acc urat ely the tail F < Fq. Going back to the 
exact result ( |4.2| ) we find for the tail of the distribution 
the expression 



P(F) 



NT^A 



[l + erf(u)] +C'(7V-3/2)^ ^4g^ 



where we have defined u = VW2(F/Fo - 1). 



The distribution Pl ( F) of the lasing mode follows from 
P(F) by means of Eq. ( |3.4|) . We find that it has a pro- 
nounced maximum at a value Mmax determined by 



exp(- 



J L-lV^{g + l) 



[l + erf(w 

n 



/2N 



(4.7) 



For L ^ y/ N (and hence also in the good cavity regime) 
we find Minax — Vln L < 0, and the deviation of F from 
Fo is of order A^/N < Fo (as long as L < e^). 



B. Mean Petermann factor 

Eigenfunction correlations of non-Hermitian operators 
have been studied in Refs. The eigenfunction 

autocorrelator considered in these studies is directly con- 
nected to the Petermann factor K. Ref. provides a 
convenient expression of the mean Petermann factor. 



MTT{K)n,rp{^) 



lim ( tr — 1- ^ 



(4.8) 



In Ref. ||l7|] this average has been calculated perturba- 
tively for 2> 1, with the result 



{K) 



(4.9) 



for Fo < F < Fo/(l - T). This result is at the same level 
of approximation as Eq. (^^) for the distribution of the 
decay rates, i. e. it does not describe the range F < Fq of 
atypically small decay rates. Since that is precisely the 
range that we need for the Petermann factor, we can not 
use the existing perturbative results. We have calculated 
the mean Petermann factor non-perturbatively for any F 
and N , assuming broken time-reversal symmetry. The 
derivation is given in Appendix The final result for 
the mean Petermann factor is 



(^>o,r = 1 



2S'(7rF/A) 



T^{nT/A)T2{TTVlAy 



(4.10a) 
(4.10b) 



with JFi and JF2 given in Eq. (^^). 

For identical gn = g 'we can use Eq. (4.3) and obtain 
by successive integrations by parts 



Siy) 



n=0 



-y 



dy-"- 



-gy 



d f sinhy 
dy\y 



(4.11) 



For iV — 1 and F ^ A we recover the single-channel 



result (3.24) of the previous section. In what follows we 
will continue to assume for simplicity that all g„'s are 
equal to a common value g. 
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FIG. 5. Average Petermann factor (iv) as a function of 
the decay rate F for different values of fully transmitted 
scatt ering channels. The solid curves are the analytical re- 
sult (i.lO), the data points are obtained by a numerical sim- 



ulation. Time-reversal symmetry is broken. 



The large- behavior can be conveniently studied from 
the expression 



S{y) 



1 



V(iv-i)!y,(,_i) 

X ~ {g - l)v\[x - {g ^ l)y\, 



(4.12) 



because the integral permits a sa ddle -point approxima- 
tion. For F > Fq we recover Eq. (4.9), but now we can 
also study the precise behavior of the mean Petermann 
factor for F < Fq, hence also for decay rates relevant for 
the lasing mode. The results will again be presented in 
terms of the rescaled parameter u = •y]V/2(r/Fo — 1)- 



We expand the integrands in Eqs. (4^) and ( 1.12 ) around 
the saddle point &t x = N (which coincides with the up- 
per integration limit at F = Fq) and keep the first non- 
Gaussian correction. This yields 



{K)nx = tV2N[F{u) +u]- T{g - 1) 

+ TF{u) [{3-g)u+ju^ + |(1 

^ exp(-^^) 
^ ' 0F[l+crf(«)]- 

For F = Fq {u — 0) this simplifies to 



u')Fiu)] 
(4.13a) 

(4.13b) 



(4.14) 



We see that the mean Petermann factor varies on the 
s am e scale of F as the decay-rate distribution P(F), Eq. 
(4.6). However, while P{T) decays exponentially for u <C 
— 1, the mean Petermann factor displays an algebraic tail 
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N 



10 



12 



FIG. 6. Average of the Petermann factor A" at F = Fq 
as function of the number N o f full y transmitted scattering 
channels. The analytical result (4.10) for broken time-reversal 
symmetry (full curve) is compared with the result of a numer- 
ical simulation (open circles for broken time-reversal symme- 
try, filled circles for preserved ti me-rev ersal symmetry). The 
dashed line is the large- Af result (4.14). The inset shows the 
distribution of A' at F = Fo for A' = 10. 



(K) 



Ty/N 
uV2 



l-T + 0{u-^). 



(4.15) 



For an amplification window ila — LA with L ^ \/ N we 
found in Section IV A| that the decay rate F of the lasing 
mode drops below Fq (the rescaled parameter Umax ~ 
— VhiL). Still, the mean Petermann factor 



(4.16) 



remains parametrically larger than unity (as long as L <C 

We now compare our analytical findings with the re- 
sults of numerical simulations. We generated a large 
number of random matrices Ti. with dimension M = 120 
{M = 200) for Af = 2, 4, 6, 8 (A^ = 10, 12) fully transmit- 
ted scattering channels {g = T = 1). Fig. ^ shows the 
mean K at given F. We find excellent agreement with 
our analytical result ( 4.10 ). 

The behavior {K) ~ ViV at F = Fq is shown in Fig. 
^. The inset depicts the distribution of AT at F = Fq for 
A^ = 10, which only can be accessed numerically. We see 
that the mean Petermann factor is somewhat larger than 
the most probable (or modal) value. 
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r/r, 







FIG. 7. Theoretical expectation (4.17) (full curve) and the 
result of a numerical simulation (data points) for the average 
Petermann factor K in the presence of time-reversal symme- 
try, as a function of the decay rate F for 2 fully transmitted 
scattering channels. 



C. Preserved time-reversal symmetry 



I 1 






N=n • 




, •10- 




. ; • ] -s" 




• I • • * , •6* 


• • • 






• * , • • • • 

• • * N = 2 " 

1 



r/To 

FIG. 8. Results of a numerical simulation of the average 
Petermann factor (A") in the presence of time- reversal symme- 
try, as a function of the decay rate F for A'^ fully transmitted 
scattering channels. 



In the derivation of the mean Petermann factor for bro- 
ken time-reversal symmetry (Appendix ^) it turned out 
that the final resuh is formally connected to the expres- 
sion for the decay-rate distribution -P(r), inasmuch as 
both expressions are built from the factors Ti (involving 
non-compact bosonic degrees of freedom of the saddle- 
point manifold) and T2 (involving compact bosonic de- 
grees of freedom of that manifold). We tried to trans- 
late this description to the case of preserved time-reversal 
symmetry (/? = 1), by operating in the same way on the 
compact and non-compact factors of the expression of 
Ref. ijl^, but could obtain a satisfactory result only for 
iV = 2, 



{K) = 



1 r(r ~ Fo) exp(F/ro) + sinh(F/ro) 



2ro r cosh(r/ro) - Fq sinh(F/Fo 



(4.17) 



In Fig. 1^ this expression is compared to the result of a 
numerical simulation. 

For larger numbers of channels we can draw our con- 
clusions from the numerical results that are presented in 
Fig. ||. Interestingly enough the data points for N chan- 
nels are close to the results for broken time-reversal sym- 
metry with N /2 channels, when the decay rate is given in 
units of Fq. This is illustrated for = 8 in Fig. Such 
a rule of thumb (motivated by the number of real degrees 
of freedom that enter the non-Hcrmitian part of 7i) was 
already known for the decay rate distribution (inset in 
Fig. H) . Hence the Petermann factor for the lasing mode 
should again display a sublinear growth with increasing 
channel number N . This expectation is indeed confirmed 




r/Fn 



FIG. 9. Average Petermann factor {K) for TV = 4, /3 = 2 
[ open circles: result of a numerical simulation, curve: Eq. 
(4.10)] and for N = [3 = 1 (filled circles: result of a numer- 
ical simulation). The parameter Fq equals NA/2n in both 
cases, so it is twice as large for (3 — 2 as for P — 1. The inset 
depicts the probability distribution of F. 
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by the numerical simulations, see the filled circles in Fig. 



i 



V. DISCUSSION 



The Petermann factor K enters the fundamental lower 
limit of the laser linewidth due to vacuum fluctuations 
and is a measure of the non-orthogonality of cavity 
modes. We related the Petermann factor to the residue 
of the scattering-matrix pole that pertains to the lasing 
mode and computed statistical properties of K in an en- 
semble of chaotic cavities. The technical complications 
that had to be overcome arise from the fact that laser 
action selects a mode which has a small decay rate F, 
and hence belongs to a pole that lies anomalously close 
to the real axis. Parametrically large Petermann factors 
oc y/N arise when the number N of scattering channels 
is large. For a single scattering channel the mean Peter- 
mann factor depends non-analytically on the transmis- 
sion probability T . 

The quantity K is also of fundamental significance in 
the general theory of scattering resonances, where it en- 
ters the width-to-height relation of resonance peaks and 
determines the scattering strength of a quasi-bound state 
with given decay rate F. If we write the scattering matrix 
(2.6) in the form 



Snni + 0-„Cr,^('^ - f2 + iT/2y 



(5.1) 



then the scattering strengths cr„, are related to F by 
a sum rule. For resonances close to the real axis (F ^ A) 
the relation is 



E 



(5.2) 



For poles deeper in the complex plane, however, the sum 
rule has to be replaced by 



E 



K > 1 



(5.3) 



The method of filter diagonalization (or harmonic inver- 
sion) that was used in Ref. to obtain for the 
molecular ion the location of poles even deep in the com- 
plex plane can also be employed to determine the corre- 
sponding residues, and hence K. 
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APPENDIX A: JOINT DISTRIBUTION OF A 
AND B 



We calculate the joint di stribu tion P{A, B) [ Eq. (P-l?!) ] 
of the quantities A [Eq. ( p^ ] and B [Eq. (|3.16|)j by 
generalizing the theory of Ref. [^. We give the lasing 
mode uji the new index M and assume that it lies at the 



center of the semicircle (2.2), lom = 0. Other choices just 
renormalize the mean modal spacing A, which we can set 
to A = 1. The quantities A and B are then of the form 



M-l 

E 



M-l I 

^ I, 



m—l 



(Al) 



The joint probability distribution of A and B, 
P{A,B) ^(sIa 



M-l I 



771 — 1 



xSl B- 



AI~1 I 

I,:' 



(A2) 



}• 



is obtained by averaging over the variables {|am| 
The quantities jamj^ are independent numbers with 
probability distribution (3.1). The joint probability dis- 
tribution of the eigenfrequencies {uj„i} of the closed cav- 
ity is the eigenvalue distribution of the Gaussian ensem- 
bles (2.1) of random- matrix theory. 



P({w,„}) cx J]^ |t^i - ujj f exp 



^ k 



(A3) 



The parameter K defined in Eq. ( 2.1g ) appears as a 
measure of mode non-orthogonality also in problems out- 
side of scattering theory. These problems inv olve non- 
Hermitian operators that are not of the form ( |2.3| ) [Q. Yi 
Many applications share the common feature that they 
can be addressed statistically by an ensemble descrip- 
tion, and that the physically relevant modes lie at the 
boundary of the complex eigenvalue spectrum. The non- 
perturbative statistical methods reported in this paper 
should prove useful in the investigation of some of these 
problems as well. 



Our choice A = 1 translates into fi = M/tt. 

The joint probability distribution of the eigenvalues 
{tOm} {m =^l...Af — 1) is found by setting cj^v/ — 
in Eq. (A_3). It factorizes into the eigenvalue dis- 
tribution of A/ — 1 dimensional Gaussian matrices H' 



again distributed according to Eq. (2.1)], and the term 
■iT'kP = |detH'|'3. 



In the first step of our calculation, we use the Fourier 
representation of the (5-functions in Eq. (A2) and write 



P{A,B) 

/ °? °? M-l 

0,1 J dx J dye"^+'^^) n J d\amfP{\a„ 



-oo — oo 
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X exp 



-IX 



Af-l 

E 

m— 1 



M-1 



(A4) 



where the average refers to the variables {wm}. The in- 
tegrals over |amP can be performed, resulting in 



det H'"^^ 

det[77'2 + 2iw{xH' + y) / tt"^ fif 



(A5) 



where the average is now over the Gaussian ensemble of 
_ff'-matrices. It is our goal to relate this average to au- 
tocorrelators of the secular polynomial of Gaussian dis- 
tributed random matrices, given in Refs. ]3^ , |3^ . 

The determinant in the denominator can be expressed 
as a Gaussian integral. 



P[A,B) 



dz / dH' det H'^f^ 



X exp 



4M 



2iw 



2 {xH' + y)]z 



(A6) 



where the M — 1 dimensional vector z is real (complex) 
for P = 1 (2). Since our original expression did only 
depend on the eigenvalues of H' , the formulation above 
is invariant under orthogonal (unitary) transformations 
of H\ and we can choose a basis in which z points into 
the direction of the last basis vector (index M ~ 1). Let 
us denote the Hamiltonian in the block form 



V h 

ht g 



(A7) 



Here F is a (M — 2) x (M — 2) matrix, g a number, and 
h a {M — 2) dimensional vector. In this notation. 



dh 



P{A, B)oc J dx J dy e'^^+'y^ J dz J dg j dV 

Xdet [y2/3(^_hty-lh)2/3] 



X exp 
X exp 



-^ig^ + 2\h\' + trV') 



2iw , . 



(A8) 



The integrals over x and y give 5- functions, 
oc Jdzjdgjdvjdh det [v'^f^{g - h^V-'^hf'^] 



X exp 



-g(5^-K2|hp + try^)-|zp(5^ + |hn 



xd{A-gB)6{B- 2w\z\^/(3tt^ 



(A9) 



We then integrate over g and z, 
P{A,B) 

A 



oc I dVdh det 

/37r2 



2/3' 



B 



|(Af-l)-2 



X exp 



4M 



(2|hp +trl/2 



Ptt'^B (A 



2w \B^ 



(AlO) 

We already anticipated B ^ 1/M and omitted in the 
exponent a term - Pn^ A"^ / AM B^ . 

The integral over h can be interpreted as an average 
over Gaussian random variables with variance 

- (l'..R^i^,^»^(i 



7r2 B/w + l/M 



MB 



(All) 

For the stochastic interpretation one also has to supply 
the normalization constants proportional to 

W ^/3(M-2)/2 

exp 



TT^B 



Pw 
'2B 



(A12) 



The integral over V is another Gaussian average, and 
thus 



P{A,B) cx QfjB^-^exp 



Pw 



Qp = ( det 



A 



1/2/5 ( - _ ht^-lh 



1-K 



2/3- 



^2^2 



(Al3a) 
(A13b) 



After averaging over h, one has now to consider for 
P=l 



}i = ( det 



.A 



V^—^ + h^V^ [(tr + 2tr y-^] 



(AM) 

where only the even terms in V have been kept. The 
ratio of coefhcients in this polynomial in A/B can be 
calculated from the autocorrelator [BSl 



{det{V + Lu){V + Lo')) 
(det V^) 
3 d siuTTX 



2x dx 



(A15) 



of the secular polynomial of Gaussian distributed real 
matrices V. This is achieved by expressing the products 
of traces and determinants through secular coefficients, 
and these then as derivatives of the secular determinant, 



(detl/2(try-i)2) 
(det V^) 

2 (det y2(trT/-2)) 



92 



duoduj 



-Gi(l^,0) 



(det F2) 



5' 
92 

4_Gi(c.,0) 



(A16a) 
(A16b) 



U1=0 
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[ We u sed t he tra nslational invariance of G{lo,lo')] Eqs. 
( IaTiI ) and (^) yield 



h « 



^2 7r2B2' 

For /3 — 2, the average over h yields the expression 



(A17) 



(Al8a) 

qi = 6{detV^[itrV-^f + tTV-'^]) , (A18b) 
+ 8try-HrV"--'^ + 6try-^ + 3(try-2)2]^ ^ ^^-,^3^,^ 

The coefRcients can now be computed from the four-point 
correlator of the Gaussian unitary ensemble |Q , 

G2(wi, ^2, W3, W4) 

(det(y + iUi){V + 0J2){V + uj3){V + 0J4)) 



3 
2^ 



(det V^) 

cos n(uJl + LU2 — LOS — UJi) 
{UJI - W3)(lJi - W4)(t^2 - ^3)(W2 - ^^4) 

COS 7r(wi + — ^2 — W4) 

(Wl - CJ2)(W1 - W4)(W3 - W2)(^^3 - ^4) 

cos7r(ti;i + W4 — 0^3 — 0)2) 



(cJl - CJ3)(t^l - t^2)(t^4 - t^3)('^4 - UJ2) 

3 

G2{(^, 0, 0, 0) = — ^(sinTrcLi — ttuj cos ttw), 
3 



G2{lu,u;,0,0) 



In this case 



-(cos27rtj - 1 + 27r w 



(Al9a) 
(A19b) 
(Al9c) 



qi = ^[6G2(w,w,0,0)-18G2(a;, 0,0,0)] 



92 = ^ [10G2(t^, t^, 0, 0) - 15G2(a>, 0, 0, 0)] 



which gives 



Q2 oc Qf. 



j=0 

(A20a) 

1^=0 

(A20b) 



(A21) 



Collecting results we obtain Eq. (3.17), where we also 
included the normalization constant. 



APPENDIX B: DERIVATION OF EQ. ( |4.10p FOR 
THE MEAN PETERMANN FACTOR 

The computa tion of the mean Petermann factor from 
the expression (4.S) is facilitated by the fact that it can 
be obtained from the same generating function [|l5|,|4| , 



_ / det[{Lu - n)iuj* - H^) - {ui - is){u2 - ie)] 
~ \ det[(tj - n)iuj* - Tit) - (u^ + i£)(u2 + ie)] 



(Bl) 



as the distribution function 
tr 



p(uj) — lim 

0+ 



X 



{uj* - H^){uj - n) + 
e 

iuj-H)iuj* -n'')+e 



(B2) 



of poles in the complex plane. (The distribution of poles 
is related to the distribution of decay rates by P(r) = 
iAp(a;)|^=o_ir/2-) The relations are 



7rp(w) 



lim 



92 



92 



52 



0+ \dLU2duj2 ^ duj2du}l 
X ^'(wi,c^2,0,0,e)^^„^^„, 
M'r:{K)^^Tp{^) 

= - hm - — — «'(w,u;,Mi,M2,e 
e^o+ 4 aui 0U2 



(B3) 



(B4) 



U\—U2—0 



Most of the analysis runs therefore in parallel with the 
calculation of p{lo) in Ref. |T^. We restrict ourselves 
to the case of broken time-reversal symmetry, where the 
algebra is less involved. 



The ratio of determinants in Eq. (Bl) can be written 
as a superdeterminant, which in turn can be expressed as 
a Gaussian integral over bosonic and fermionic variables. 



^'(t^,w,e,ui,U2) = {-lf'{Sdet-\A)) 
= (-1 



/ I d^^ I d*e'*'^* ) . (B5) 



The matrix A is 



A 



/ uj-H ie + ui \ 

uj — H —ie + ui 

-ie-U2 -UJ* +T-0 

\ -ie + U2 UJ* -rO J 



{n-H)®L + i [t^W'^W - - ) (g) fJ^L 
— ie ® UxL + u&xL. 



(B6) 



The vector * = *! ® \I>2 © *3 © *4 is a 4M- 
dimensional supervector consisting of two M-dimensional 
bosonic entries '4'q. with a = 1 and 3, supplemented 
by two Af-dimensional fermionic entries with a — 2 
and 4. We encounter the 4-dimensional supermatrices 
L = diag(l, 1, — 1, 1), u — diag(— iti, ui, — it2, U2), and 
(7j — (ji 1^ IL2, where cr^ are the usual Pauli matrices [e.g. 
=diag(l, 1,-1,-1)]. 

The linear appearance of H in the exponent of Eq. 



(B5) facilitates the ensemble average with the distribu- 
tion function (|]^), since the integral over the indepen- 
dent components of H factorizes, and each single integral 
is Gaussian. The result is 
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exp[-i*^i? (g) L*] ) = exp 



M 



(B7a) 
(B7b) 



The order of R in the exponent is reduced from quadratic 
to hnear by a Hubbard- Stratonovich transformation, 
based on the identity 



exp 



dS exp 



Str(Li?)^ 



(B8) 



The integral over and VI'^ is again Gaussian and results 
in 



dS exp 



-A/Str ( y + lnS' 



Sdet-i (1 + C) 



(B9a) 



(n + i (ttW^W 



I az — iecTx + udx 



(B9b) 

One now can write Sdet"^(l+C) = exp[-Strln(l+C)] 
and expand the logarithm to first order in F, e, and the 
source term J; in addition we set 17 = and pass from 
the generating fu ncti on to the mean Petermann factor 
according to Eq. (B4). This gives 



MTr(K)n.rp{uj) 



- MStr 




■ y 



— StrCT^^* ^ + I— Str(Ta;5' ^ 



N 

S-^ [| Sdet^i (I4 



(BIO) 



The traces iiijA = An + Ajj operate only on the indi- 
cated subspaces. We introduced the rescaled variables 
y = -27rIma;/A = ttT/A and e' = 27r£/A. In what 
follows we will write e instead of e'. 

The condition Af 3> 1 justifies a saddle-point approxi- 
mation. The main contribution to the preceding integral 
comes from points for which the first part of the exponent 
is minimal, that is from the solutions of 



4 + S' = 
S 



= -1. 



(Bll) 



With S ~ iQ, the solutions fulfil l = 1. As inherited 
from the definition of R in Eq. ( B7b| ), QL is a Hermi- 
tian matrix and Q — T^^Q^iagT can be diagonalized by 



a pseudounitary supermatrix f G f/(l,l/2) (these ma- 
trices fulfill T'^LT = L). The largest manifold which 
respects the definiteness-requirements on Q is obtained 
by the choice Qdiag = (^z- However, rotations in the 
block a = 1,3 and in the block a = 2,4 leave Q in- 
variant; the saddle-point manifold is hence covered ex- 
actly once if we take the T matrices from the coset space 
(7(l,l/2)/C/(l/l)xt/(l/l). 

A convenient parameterization of the coset space has 
been given by Efetov [Ba , 



T = 











exo, idiag(0i,z02) 

''''PUdiag(0i,z02) 



U 
V 



-pp*/2 
P* 1 



p*p/2 



V = 



ia* 1 



(tV/2 



(B12a) 



(B12b) 



(B12c) 



with bosonic variables ^i, ^2, 0i, and 02, and fermionic 
variables p, p*, cr, and a* . We introduce Ai = cosh 6*1 
and A2 = cos 02- In this parameterization 



(B13a) 

- ipa*\ 



Strang = 2(Ai-A2), 

Stro-^Q = -sinh6'ie*"*i[(l + pp* I2){\ - aa* 12) 
^ sinh6'ie-''^4(l + pp* /2){1 - aa* 12) - iap*] 
+ i sm02e"l'^ [(1 + p*p /2){1 - a*a/2) - ip*a] 
~ i sineae-"^^ [(1 + pV/2)(l - a*a/2) - ia* p], (B13b) 
tri2crx(3 = - sinh6lie*"*i [(1 + pp* /2){1 - aa* /2) 
- i sin 6*26^'^' [(1 -I- P*p/2){1 - a* a 12) - iap*], 
tm<7xQ = sinh^ie-*"^! [(1 + pp* /2){1 - aa*/2) - 



- i<7*p] 
(B13c) 



ip a\ 



+ i sin 6*26-''^^ [(1 -I- p*p/2){l - a*a/2) - ipa*] , (B13d) 



Sdet~^[l4 + WnazQ] 



Qn + A2 



9,1 + Ai 
The integration measure is 



dQ = 



dXi d\2 d(j)i d(j)2 dp* dpdia* dia 
(27r)2(Ai-A2)2 • 



(B13e) 



(B14) 



In order to integrate over the fermionic variables we 
have to expand in these quantities and only keep the term 
in which all four variables appear linearly. The angle 02 
appears in the pre-exponential factor as well as in the 
exponential term exp(— esin02 sin02)- We expand the 
exponential and integrate over 02- Only terms of order 
e" sinh™ 9i with n < m survive the limit £ ^ 0. We 
discard all other terms and obtain 

-4^(X)o.rp(c^) 
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lim 



X exp 



27r 



1 sin(/)i + ?/(Ai - A2) 



d<h_ 

N 

n 



gn + A2 



(B15a) 



D = -ie sinh di sin (2 sin^ 6*2 + 1 sinh^ 6^1 ) 



— sinh^ 9i [sinh 



^ ^1 cos^ (pi 



(3 cos^ 



5 sin^ sin^ ^2] 



is sinh 0i sin 0i sin 82 ( 
j^e*^ sinh^ 9i sin^ 0i sin^ 02 ■ 



S- sin^ </)i 



(B15b) 



It is convenient to bring the factor D into a form which 
involves only in the combination zi = ~i sinh 0i sin 0i , 
because such terms can be expressed as derivatives with 
r espec t to e of the exponential exp(£Zi) appearing in Eq. 
( B15a ). This goal can be achieved by integrating by parts 
all terms that involve cos 01. Effectively this amounts 
to the substitutions esinh^i sin(/>i cos^ ipi i(sin^ — 
cos^ (f>i) and e sinh^i cos^ 0i — s- i sin0i, resulting in 



D = ezi(| sin^ 6*2 + 2sinh^ 6*1) + ^e^zf sin^ ( 
- ie^z? sin^ 02- 



(B16) 



Mathematically these expressions are quite similar to 
those obtained for the decay-rate distribution in Ref. . 
By a simple substitution rule that relates to each other 
the terms of different order in e, we now rewrite D in 
a way that allows to make direct contact to Ref. 
yieldin g a r esult in terms of the two functions J^i^2 given 
in Eq. ( |4.2| ): As in Ref. we express the factors {gn + 
Ai)~^ as an integral of exponential functions 



1 

— — = / dsn exp[-s„(.g„ 
+ Ai Jo 



Ai)]. 



(B17) 



We also write (Ai — A2)^ = dx x exp[—x{Xi — A2)]. 
Then the integrations over 0i and 4>i can be performed, 
and e only appears in a factor 



exp[ 



y 



(BI8) 



with y' ^ y — X — s„. The limiting value for e — s- of 
the derivatives 

£"^$(£, y') ^ CnSiy'), Ci = -C2 - C3/2 = -2, 

(B19) 

2ezi 

and e^zf — £zi, which gives D = 2ezi{Xi — A|). As a 
result, we obtain 



amounts in Eq. ( |B16| ) to the substitutions e 



3^3 



A 



P{T) 



{K)nM = /o(vrr/A) + 2/i(7rr/A), (B20) 



1 d f°° 
/.(y) = --lhne-y^ dX,J_^ 



dXn 



(Ai - A2)- 



N 



Jo(ejA2-l)exp[2;(Ai-A2)]n 



gn + A2 
\ gn + Al ' 



(B21) 



where Jq is a Bessel function. By comparing expressions 
with Ref. we recognize that Ia{y) — !Fi{y)!F2[y) = 
(A/7r)P(r = A^/tt) [cf. Eq. (jij)], while 



h{y) 



j\y' Ti{y')^T2{y'). (B22) 



This concludes the derivation of the final result (4.10). 
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